#include "TMath.h"
class Event{
    public:
        Event(){ }
        Event(
                vector<vector<double> > &_theta,
                vector<vector<double> > &_phi
             ){

            vector<double> cosvalue, sinvalue, value;

            costhetas.clear(); sinthetas.clear(); thetas.clear();
            for(int i=0; i<_theta.size(); i++){
                cosvalue.clear(); sinvalue.clear(); value.clear();
                for(int j=0; j<_theta[i].size(); j++){
                    cosvalue.push_back(TMath::Cos(_theta[i][j]));
                    sinvalue.push_back(TMath::Sin(_theta[i][j]));
                    value.push_back(_theta[i][j]);
                }
                costhetas.push_back(cosvalue);
                sinthetas.push_back(sinvalue);
                thetas.push_back(value);
            }

            cosphis.clear(); sinphis.clear();phis.clear();
            for(int i=0; i<_phi.size(); i++){
                cosvalue.clear(); sinvalue.clear(); value.clear();
                for(int j=0; j<_phi[i].size(); j++){
                    cosvalue.push_back(TMath::Cos(_phi[i][j]));
                    sinvalue.push_back(TMath::Sin(_phi[i][j]));
                    value.push_back(_phi[i][j]);
                }
                cosphis.push_back(cosvalue);
                sinphis.push_back(sinvalue);
                phis.push_back(value);
            }
            numEvt_total = _phi.size();
        }


        vector<vector<double> > thetas, phis;
        vector<vector<double> > costhetas, sinthetas;
        vector<vector<double> > cosphis, sinphis;

        double numEvt_total;
};
